Energy relaxation in the Gross-Pitaevskii equation 
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We introduce a dissipation term in the Gross-Pitaevskii equation that describes the stimulated 
relaxation of condensed bosons due to scattering with a different type of particles. This situation 
applies to Bose-Einstein condensates of quasi-particles in the solid state, such as magnons and 
excitons. Our model is compatible with the phenomenology of superfluidity: supercurrents are 
stable up to a critical speed and decay when they are faster. 
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The Gross-Pitaevskii equation (GPE) provides an ex- 
cellent description of an isolated dilute Bose gas at low 
temperature. For example, the GPE accurately describes 
the density profiles and frequencies of elementary excita- 
tions of trapped ultracold atomic bose gases [lj. 

Systems featuring Bose-Einstein Condensation (BEC) 
of quasi-particles in the solid state, such as (exciton-) 
polaritons [2] and magnons [3[ are however not so well 
isolated from their environment as the atomic gases. In 
semiconductor microcavities for example, the polaritons 
that form a Bose Einstein condensate undergo scatter- 
ing with lattice phonons and with high energy excitons. 
These interactions allow the polaritons to relax toward 
lower energy. 

The Boltzmann equation is the easiest way to model 
these scattering dynamics theoretically [4]. The semi- 
classical approximation involved in this approach be- 
comes however problematic when the polariton gas enters 
the condensed phase. Most importantly, supercurrents, 
that are a gradient of the phase of the polariton field, are 
missed. This precludes e.g. the description of vortices 
in polariton condensates [5] , that appear naturally in the 
GPE. 

In order to allow for a simultaneous description super- 
currents and of relaxation due to interactions with the 
environment, we will derive in this Letter a dissipation 
term that can be added to the GPE. This term models 
transitions between components of the GPE field at dif- 
ferent frequencies (see Fig.[T]). For specificity, we will use 
in the following the terminology of 'polaritons' relaxing 
through 'phonon scattering', but our formalism is appli- 
cable to any dilute Bose gas that dissipates energy into 
the environment. 

An important requirement of the theory is compati- 
bility with superfluidity. According to the Landau crite- 
rion [1], the superflow cannot relax due to phonon scat- 
tering as long as the superfluid velocity is below the criti- 
cal speed. A Boltzmann description that includes phonon 
scattering violates this requirement: a Bose gas with ini- 
tially all the particles at finite momentum k = k c will re- 
lax to a state with the majority of particles in the ground 
state k = 0. We will show below that our GPE based 
model for phonon relaxation is instead fully in agree- 
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FIG. 1. Schematic illustration of the relaxation of polaritons 
due to scattering with phonons. Several modes in a random 
potential are macroscopically occupied. Where they spatially 
overlap, relaxation takes locally place from the high to the 
low frequency modes, as indicated by the arrows. 



ment with the Landau criterion for stability of superflow 
slower than the critical speed. In addition, it provides 
a dynamical model for the dissipation of superflow when 
the speed exceeds the critical one. 

A kinetic model formulated in terms of transitions be- 
tween states characterized by their frequencies involves 
a coarse graining of the time, analogous to the coarse 
graining of space used in the semiclassical Boltzmann 
equation [6[ . We will use the following definition of the 
time dependent spectrum of the bose field ip(x,t): 
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where T is the coarse graining time step. The inverse 
transform reads 
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Spontaneous scattering is a quantum feature that is dif- 
ficult to describe with a classical field model such as the 
GPE. We therefore restrict to the stimulated relaxation 
processes: 



dn UJ (x,t) = n u) (xjt) V^ r(cj,cj / )n u; /(x, t) dt, (3) 



where r(uo,uo f ) is the net scattering rate from the mode 
at frequency uo' to the mode at a;. In order to conserve 
particle number, the relaxation rate should change sign 
under exchange of uo and uo': r(uo,uo') = —r(uo',uo). The 
first term in the expansion of the rate as a function of 
the frequency difference is therefore 



r(uo,uo') = n (uo' — a;), 



(4) 



where the relaxation constant k has the dimension of an 
inverse density. Eq. (J2J) is for example the form of the 
scattering rate obtained with the golden rule for the re- 
laxation of quantum well excitons due to scattering with 
acoustic phonons [7j. 

Under the change of the density (|3]) , the wave function 
varies as 
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(5) 
where we used the fact that a stimulated relaxation pro- 
cess does not change the phase of ip. Expanding Eq. (|5j) 
to first order in dn^, substituting Eqns. (J3j) and (|4]) and 
using the inverse transform (|2]), one obtains for the re- 
laxation dynamics of the wave function 



dij)(x,t) K,n(x,i) 
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In the derivation, boundary terms in the partial integra- 
tion that scale as 1/T were neglected, because we con- 
sider the limit of a large coarse graining time where 1/T 
is smaller than any other relevant frequency scale. In 
Eq. (j6]) , n and j2 are the time averaged density and chem- 
ical potential respectively: 
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The real part is taken in Eq. (j8j), because it is readily 
shown with the Madelung transformation ip = y/ne 10 , 
that the imaginary part of the integral in Eq. (J5J) scales 
as 1/T. Eq. ((6]) gives the term that was sought to de- 
scribe the stimulated relaxation of a classical Bose field 
due to local scattering with phonon like particles and is 
the central result of this Letter. 

The right hand side in Eq. (J6j) resembles the frequency 
dependent amplification term that was introduced in 
Ref. [8] to describe an energy dependent gain mechanism: 
dip/dt = (P/VLk)[Qk — id/dt]ifj. Two main differences 
should be emphasized. First, here the relaxation term 
is proportional to the polariton density n, whereas in 
Ref. [8], the relaxation is proportional to the gain from 
the reservoir P. The second important difference is that 
the gain cutoff frequency Qk is replaced by the average 



frequency ft. As a consequence, particle loss and gain 
balance each other and there is no net gain in Eq. (J6j). 
The modes with frequency uo > fl experience loss, where 
the ones with uo < ft are amplified. Phonon absorption is 
thus not included in our model. The relaxation (J6]) is due 
to the interaction with an environment that is effectively 
at zero temperature. 

As a first application, we add the dissipation term (|6j) 
to the GPE for noninteracting particles in two coupled 
levels. In the context of polariton condensation, the cou- 
pled two-level system describes for example the situation 
of spatially distinct potential wells [9] or the polarization 
degree [10]. The full dynamics is 
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for j — 1,2, where J is the coupling energy. Fig. [2] shows 
the time evolution of the densities in the two states, ob- 
tained from the numerical integration of Eq. ((9]). The 
initial wave function is taken to have equal amplitudes 
on the two states and a small phase difference between 
them. Due to the dissipation, the Rabi oscillations are 
damped. Note that the time evolution (J9j) drives the 
system automatically into the ground state and that the 
eigenstates of the Hamiltonian were not explicitly needed 
to compute the phonon scattering rates. The contribu- 
tion of the excited state to the wave function is shown 
in the inset of Fig. [2j Its exponential decay is excellently 
reproduced with the rate obtained from Eq. (J3]). 

We have numerically checked that the dissipative dy- 
namics conserves particle number in the limit of a small 
time step in the numerical integration. We correct the 
residual deviation by restoring the norm of the wave func- 
tion after each application of the dissipation operator. 

We now show that in the case of the interacting ho- 
mogeneous Bose gas, our description of the dissipation 
is fully compatible with the phenomenology of persistent 
supercurrents. The Gross-Pitaevskii equation for the ho- 
mogeneous interacting Bose gas supplemented with the 
phonon scattering term reads 
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The usual steady state solution ipo(x, t) = ^/n^e~' lflt +' lkcX 
is still a solution of Eq. ([TO]) with ji = k^/2m + gn Cl 
where k c is the condensate momentum, g the interaction 
constant and n c the condensate density. 

The excitations on top of the condensate can be de- 
scribed by the wave function ip(x,t) = i/jo(x,t)[l + 
ue ikx-iu(k)t _|_ v * e -ikx+iuj* {k)t^ Linearizing the equations 

of motion in u and v is not more difficult than in the 
standard case, because to first order in u and v, the av- 
erage density n and frequency j2 do not change. In the 
small dissipation limit ktl <C 1, the frequencies of the 




FIG. 2. Dynamics of the relative phase (A0) and density 
on the first well (m) for two wells coupled by tunneling and 
dissipating through the emission of phonons. The initial con- 
dition was taken to be ^1,2 = e * 0-05 Inset: The decay of 
the projection of the wave function on the excited state is ex- 
cellently reproduced by the exponential decay exp(— 2 J Knit) 
(dashed line), where the phonon relaxation constant was 
taken k, = 0.05 and the coarse graining time scale T = 20 J -1 . 



elementary excitations read : 



uj(k) = ±^uj%(k) - (ftn/2) 2 + kk c /m 
k 2 /2m + gn 
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[uj B (k)±kk c ]. (11) 



Note that the correction quadratic in nn to the usual Bo- 
goliubov dispersion uJB(k) = \J gnk 2 /m + k 4 /4m under 
the square root on the first line in Eq. ([TT]) is necessary to 
satisfy the Goldstone theorem that requires at least one 
branch of elementary excitations to have exactly zero real 
and imaginary part of the frequency for k —> 0. 

Importantly, for a condensate wave vector below the 
critical one fc cr i t = min[cjB(fc)/fc], the elementary excita- 
tions have a negative imaginary part and the condensate 
is dynamically stable against decay into a lower energy 
state. As soon as k c > fc C rit> the imaginary part of the ex- 
citations becomes positive and the condensate becomes 
unstable with respect to decay toward a lower momen- 
tum. In agreement with the thermodynamical consider- 
ations made by Landau [1], the phonon scattering is able 
to dissipate the supercurrent only when the condensate 
velocity is above the critical speed. In the present dy- 
namical model ([TO]) , the energetic instability gives rise to 
a dynamical one, because the dissipative environment is 
explicitly included and describes the kinetics of the su- 



percurrent's decay. 

As a last illustration, we describe the relaxation of 
polaritons in a harmonic trap as observed in Ref. |ll|. 



When polaritons were created off center, relaxation to- 
ward the lower energy state in the middle of the trap was 
observed. This experimental configuration is illustrated 
in Fig. [3] Another striking evidence of the relaxation 
mechanism modeled here is the very recent experimental 
result reported by Wertz et at. [12j. A detailed analysis 
of the latter experiments is under way. 

In polariton condensates, a steady state can be reached 
when new particles are continuously injected by optical 
excitation. In order to avoid the pinning of the conden- 
sate phase by the excitation laser, the additional particles 
are created at a frequency different from the condensate 
one (nonresonant excitation). In the case of a far de- 
tuned laser, the excitation creates free electron hole pairs 
that relax to form an excitonic reservoir. The further re- 
laxation of the excitons into the lower polariton branch 
then provides a gain for the condensate. A steady state 
is reached when gain and losses compensate each other. 
The stimulated part of the scattering into the lower po- 
lariton branch can be modeled by introducing a gain term 
in the Gross-Pitaevskii equation 13|,LL4|. The full dynam- 
ics of the polariton gas including trapping, nonresonant 
pumping and relaxation is given by 
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Here 7 is the loss rate and the function R(n R ) describes 
the gain of the lower polariton branch thanks to stimu- 
lated scattering from the nonresonantly excited excitons 
in the reservoir with density n R . The exciton density can 
be described by the rate equation 

j t n R = - lR n R - R(n R M 2 + P, (13) 

where 7^ is the reservoir relaxation rate and P the non- 
resonant pumping rate. The effective potential V e s = 
Hrap + V exc consists of the trapping potential Vt rap and 
the blue shift due to the high energy excitons, propor- 
tional to the pumping term V exc = QP [15| . 

It is worth pointing out that Eq. ([T2]) actually describes 
two relaxation mechanisms: (i) from the excitonic reser- 
voir into the lower polariton branch and (ii) from high to 
low energy polariton states. Only the second relaxation 
process is modeled with the new dissipation mechanism 
derived here [term on the second line of Eq. ([12]) ]. 

In the absence of the relaxation term, the 
model ([T2]) , ([13]) predicts a condensate at a single 
frequency [see Fig. [3^a)]. The frequency of the conden- 
sate coincides with the potential energy at the center of 





tions. 

Physically, it is likely that a large contribution to the 
relaxation within the lower polariton branch comes from 
the scattering of polar itons with reservoir excitons. In- 
deed, estimates based on the golden rule for exciton- 
phonon scattering [7[ yield a relaxation constant of the 
order of n ~ 10 _4 /im 2 . Our numerical simulations 
show relaxation toward the bottom of the lower polari- 
ton branch for nn > 0.1. If exciton-phonon scattering 
was the only relaxation mechanism, it would require a 
polariton density n = 10 3 /im -2 , much larger than the 
experimental estimates [11] . 

Further applications of the relaxation in the descrip- 
tion of polariton condensates include the relaxation of 
nonresonantly excited polariton condensates in periodic 
potentials [16|,LL7j, in one dimensional polariton wires [12| 
and the interplay of phonon relaxation with parametric 
scattering [18[ . Our model may also be useful to describe 
the interaction of the thermal cloud with the condensed 
part of an equilibrium Bose gas at finite temperature [19| . 

We are indebted to B. Pietka, I. Carusotto, K. 
Lagoudakis, T. Liew and F. Manni for stimulating dis- 
cussions. 



FIG. 3. Relaxation of a polariton condensate that is nonres- 
onantly created by a pump profile (dashed line) that is not 
centered at the minimum of the effective trapping potential 
(full line). Panel (a): In the absence of a relaxation mech- 
anism (k = 0), the polaritons accelerate ballistically when 
moving toward the bottom of the potential: The energy re- 
solved real space distribution shows a single frequency for the 
condensate. Panel (b): When the dissipation mechanism is 
present (k, = 0.05/im 2 ), the polaritons relax to the bottom 
of the trap. The condensate frequency in the center of the 
trap is different from the condensate frequency in the pump- 
ing area. The red dashed line shows the excitation intensity 
P and the full black line represents the effective potential V e s. 
Other parameters: fay — 0.05meV, rriLp/h 2 = meV -1 /im -2 , 
HP — lmeV/xm -2 , Hjr = 0.5meV, HR(nR) — (lmeV/xm 2 )nR. 



the pumping spot. Away from the pumping spot, the 
potential energy is converted into kinetic energy. 

On the other hand, when the relaxation is included 
[Fig 09(b)], a condensate also appears at an energy close 
to the bottom of the trap. It is important to highlight 
that the steady state solution consists of two condensates 
at well defined frequencies. As our mean field model 
only describes the coherent part of the bosonic field, the 
condensate line width vanishes. 

Fig. 2 (b) is globally in agreement with the experimen- 
tal observations in Ref. (ll|. Differently from our sim- 
ulations though, in experiments polariton luminescence 
was observed at all frequencies between the pump region 
and ground state [U| . Possibly microcavity disorder and 
spontaneous relaxation are responsible for these devia- 
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